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Abstract. This letter aims at justifying the stochastic equations in terms of 
the number density variable, which are still controversial, via complementing 
Dean's approach [Dean D S 1996 J. Phys. A 29 L613]. Our course is twofold: 
First, we demonstrate that standard manipulations straightforwardly transform 
the stochastic equation of density operator, derived by Dean, to the Fokker-Planck 
equation for the (c-number) density distribution functional P{{p},t). Moreover, 
we verify the associated static solution of P{{p}, t) with the help of the conditional 
grand canonical partition function. 



1. Introduction 

In supercooled liquids, due to the dense packing and strong correlation of the 
constituent particles, the nonvibrational diffusive motions take much more time than 
collisions. In other words, the momentum and the energy flow much more quickly 
via collisions through the system than the slowly decaying number density. Recently 
there has been a considerable effort to describe such slow dynamics in liquids . We 
shall in particular discuss one of the approaches, the following stochastic equations in 
terms of the number density field p{x,t) [2-5]: 



+ ax,t), (1) 

p{x,t) 



or its equivalent, i.e. the Fokker-Planck equation for the probability distribution 
functional P{{p},t): 

In equations and (j^), the Hamiltonian H is of the free energy functional form as 

H{{p})^^ J dxdyp{x)V{x-y)p{y)+T J dx p(x)\og p{x), (3) 

L[p{x)] is the kinetic coefficient written as L[p{x)] ~ p{x,t)T with T being the 
mobility of particles, and ^ is the divergence of a random force and its correlation 
function is given by 

iax, miy, t')) = 2TV, • L[p{x) ]^yS{x - y) 5{t - t'). (4) 

For the explicit representation of this averaging, see equation (0) below. 
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These equations have attractive features for studying the slow dynamics: One is 
the physically clear incorporation of the thermally activated hopping processes via the 
last random term ^ on the right hand side (rhs) of equation (^. Also of interest is 
the density dependence of the kinetic coefficient, L oc p, that produces the nonlinear 
term of dynamic origin |0 and implies the relevance of these equations to describing 
dynamical heterogeneity . 

Nevertheless the stochasticity for the number density variable is still controversial. 
To see this, let us first mention one of some attempts [2-4], presented by Dean to 
justify the above stochastic equations (|l|) and (^: Consider here a canonical system 
of N particles interacting via a pairwise potential V{x) and surrounded by a thermal 
white noise heat bath. The ith particle then obeys the Langevin equation, 

dX (f) ^ 

= -r J2v,v[x,{t)~x,{t)] + v^it), (5) 

where the components of the noise rjf are taken to be uncorrelated as {Vi {t)rif {t}) = 
2Tr6ijS°''^6{t — t'). Dean shows, using the Ito prescription for the change of variables 
and summing over the i, that equation (|^) is transformed to the equation of the density 
operator, /3(x,t) = J2i Pi{x,t) = J2i ^[^ - ^iit)]- 
dp{x,t) 



= V • p{x, t) I dy p(y, t)VV{x - y) + TV^p{x, t) + ^ (x, t) 



iix,t) = - ^V- [p,(x,t)77,(i)]. 



(6) 



Since one finds 



dy {p{x, t)p{y, t))VV{x ~ y) + TV^Pavi^, t), (9) 



t)iiy, t)} EE I dv, iix, t)iiy, t) exp {- j dt (7) 

= 2TV,-L[p{x)]VJ{t-t'), (8) 

equation ([^) is verified so long as the operator p reads p of c-number. 

As expected, though, Marconi and Tarazona (MT) j|] subsequently object to the 
last supposition: they claim that p is to be defined by averaging p over the noise as 
pav = (p), where the subscript av is appended for emphasizing the procedure. As a 
consequence dynamical density functional equation becomes deterministic: 

dpay{x,t) ^ 

dt 

whereby the Boltzmann distribution of number density is assured as time-independent 
solution |,|9). 

To settle such controversy over the stochastic density functional equations (|l]) 
and this letter aims at complementing Dean's argument from (||) to (^) so that 
the above criticism by MT may become invalid. Our strategy is twofold: First, in the 
next section, wc demonstrate that standard manipulations transform equation 

of the density operator to the Fokker-Planck equation Moreover we verify 
in section 3, with the help of the conditional grand canonical partition function, the 
static solution Pniip}) of (§): 

Fo(M)«exp(-/3i/), (10) 

where (3 — T^^ . In the final section, to clarify the connection between the stochastic 
and the deterministic equation, we confirm using the WKB-like approximation 
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to the Fokker-Planck equation (|^) that the noise-averaged deterministic equation 
corresponds to that for the saddle-point path of P{{p},t); this reveals that MT's 
argument produces only the mean-field equation and not the first member of the 
BBGKY hierarchy | p2[ including two-point equal-time correlation function. 

2. From equation to the Fokker-Planck equation (^) 

Turning our attention to the functional space, we immediately find that the density 
operator p may be directly mapped to the distribution functional P({p},<) as 

P{{p},t) = /^llS[p{x,t) - pix)]^ (fl) 

not via the averaging pav — (p); essentially only this definition has dissolved the MT's 
critique. 

Let us then exhibit below that P{{p},t) with equation ^) satisfies the Fokker- 
Planck equation (||). We first differentiate (|ll| ) with respect to time, and obtain 

dP{{p},t) _ [ , I dp(x,t) 5 



= /rf^ ( ^^g/ J ^) mx, t) - p{x) ] n my, t) - piv) ] ) 

V-p(x,t) / dy p{y,t)VV{x - y) + TV^p{x,t) + i{x,t) 



dx 



X 



—6[p{x,t)-pix)]Y[S[piy,t)-piyy_ 



6p{x,., 



We may replace p by p using the (5-function, and hence equation ( p^ reads 
dP{{p},t) 



dt J 5p{x 



— — I dx 



V • p{x) J dy piy)VVix ~y) + T V'p{x) 
■ J dx I i{x, t) 5[p{x, t) - p{x) ] 



P{{p],t) 

(13) 



With use of the identity for an arbitrary function f ({??"}), {F{{'nf})'nf{i)) = 
2TT {5F{{r}f})/5r]f{t) ), the last bracket term on the rhs of ( [l^ ) is further transformed 
to 



V 



= Tr^^V,.V,[p,(x,i)p.(y,t)]^^^n<5[p-p]^, (14) 

where the superscript a denotes the component of x and rji, and use has been made 
of 

Spi{y,t) 1 dp.,{y,t) 

Srjf ~" 2 dyf ^ ' 

that is obtained from standard mathematical manipulation of the discretized Langevin 
equation |l^. Also, noting that relation pi{x,t)pi{y,t) = 5{x — y)pi{x,t) gives 

Va; ■ Vy[pi{x,t)pi{y,t)] = -Vx ■ Pi{x,t)Vx5{x - y) (16) 
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and replacing p by p as before, the bracket term finally reads 



iix, t) Yl 6[pix, t) - p{x) ])=TT^-p{x,t)V 



5P{{p}.t) 



(17) 



5p{x,t) 

Equation ( p^ with this is none other than the Fokker-Planck equation (||). 

Thus it has been demonstrated that the stochastic equation of density operator 
(^ leads to the Fokker-Planck equation (||) for the (c-number) density distribution 
functional [or equation (1^)]. 



3. Verification of the equilibrium distribution functional (10) 

We arrive at a time-independent solution of the Fokker-Planck equation in the large 
time limit: Po{{p}) — ^^"cat— >oo P{{p} it). Therefore it is plausible to suppose that 
the noise averaging in calculating Pq becomes equivalent to the configurational one in 
equilibrium: 

N 

^o(M)cx— n / dX,\{5[p{^,t)-p{x)] exp[-/3^F(X,-X,)]. (18) 

i—l X i^j 

The problem is then how to derive expression ( p^ ) from the above conditional partition 
function. 

Let us once move to the grand canonical system where we are to consider 



(19) 



JV=0 



with = In A being chemical potential. Introducing the auxiliary field ■0 as 
S[p{x) — p] ~ J dip exp [ iipip — p)], the configurational representation of Pq given by 
(^ and (|l|) reads 



Dtp exp 



-(3 I dxdy^ p{x) V{x - y) p{y) + ip{x) ijj{x) - e*'''(^)+^ 



(20) 



where Dtp is formally defined as Y\^dip{x). Since there is no contribution to of 
principal quadratic fluctuation of the auxiliary field ip around the saddle point path 
ipsp as shown elsewhere fl^ j, Gaussian approximation for tp reduces the functional 
integral form (poh to 



P^ (X exp 



l3H + P / dx p{x) + fip{x) 



(21) 



as found from substituting p — e'^^^p+i^ into (pO|). 

For returning to the canonical system, we have only to perform the following 
contour integral, 



Pn - 



1 



2^z/-"A^+i' ^^^^ 
where A is now a complex variable. This relation with use of the Cauchy's integral 
theorem gives back the canonical form: 

~ I f .. 1 



Pq oc exp 



-PH + (3 J dxp{x) ^^j^dX 

^-I3{H-N) if Jdxp{x) = N 

otherwise. 



^l + [N-Jdxp(x)] 



(23) 
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The static solution (^0|) has been thus verified, and supplementing Dean's discussion 
has been completed. 

4. Discussion: connection with the deterministic equation (|^) 

Here we would like to first confirm the saddle-point path of the Fokkcr-Planck (||) 
from exploiting the WKB-like approach made in Ref. Setting, similarly to the 

WKB approximation, that 

P({p},t)«exp[-/3$({p},t)], (24) 

we obtain the Hamilton-Jacobi like equation: 



a$(M,t) _ f ,6<i>i{p},t), 



Sn{p},t) SH{{p}) 



5p{x) ^p{x) 



(25) 



A short-cut way of deriving from this the most probable (or saddle-point) path {p} is 
to expand $ and H around {p} as 

$({p}, t) - $({p}, t) + ^ fdxdy [pix) ~ -p{x) ] $" [x ~ y) [p{y) ~ p{y) ] + ■■■ 

(26) 

[p{x) - p{x)]^ , 

{p}={p} 

with J{{p}) = L[p{x)y\/6H{{p})/Sp{x). Substitution of these into equation ( p5| ) 
yields in 0[p{x) — p{x) ] 



Jiip}) = Jiip}) ' 



6p 



= V • L[p[x,t)\V' 



at wj 



(27) 

{p}={p} 



Since {p} is to be in accord with the noise averaged density {pav}, equation ( p7| ) 
implies that {p{x,t)p{y,t)) = Pav{x,t)paviy,t) for the first term on the rhs of (^, 
i.e., no spatial correlation of noise-averaging. In other words, the above derivation 
reveals that the noise-averaged equation (^) is not the first member of the dynamical 
BBGKY hierarchy unlike the proposal by MT, but is only the mean-field equation for 
the saddle-point path of P{{p},t). 

We have thus validated the stochastic density functional equations, which must 
be a powerful tool for the understanding of supercooled fluids and glasses, via 
proving the irrelevance of MT's objection to Dean's argument in three ways: (i) 
demonstrating that standard manipulations enable to replace with the c-number 
density field p the corresponding operator variable p in the stochastic equation (^ 
derived by Dean, (ii) verifying the static solution (^0|) of the Fokker-Planck equation 
for the density distribution functional with the help of the conditional grand canonical 
partition function, and (iii) pointing out that the noise averaged path satisfying the 
deterministic equation (|^) merely corresponds to the saddle-point one. 

The next problem is how to solve these dynamical density functional equations, 
stochastic or deterministic. In previous works ||, ^ the static density functional 
theory has been exploited as input, and some justifications have been also described 
by Kawasaki and MT |||, js). However, the present discussion does not support these; 
from our point of view, what to suppose for incorporating the static theory remains 
an open problem. 
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